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Abstract. We study the long-time, large scale transport in a three-parameter family of 

isotropic, incompressible velocity fields with power-law spectra. Scaling law for transport 

is characterized by the scaling exponent q and the Hurst exponent H, as functions of the 

parameters. The parameter space is divided into regimes of scaling laws of different functional 

O ■ forms of the scaling exponent and the Hurst exponent. We present the full three-dimensional 

t^ . phase diagram. 

r-| '■ The limiting process is one of three kinds: Brownian motion {H = 1/2), persistent 

Oh| fractional Brownian motions (1/2 < H < 1) and regular (or smooth) motion {H = 1). 

We discover that a critical wave number divides the infrared cutoffs into three categories, 

critical, subcritical and supercritical; they give rise to different scaling laws and phase di- 

<yQ • agrams. We introduce the notions of sampling drift and eddy diffusivity, and formulate 

. variational principles to estimate the eddy diffusivity. We show that fractional Brownian 

\^ I motions result from a dominant sampling drift. 
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1 Introduction 

The movement of a passive scalar in a turbulent flow is described by the stochastic differential 
equation 

dx{t) = V{x{t),t, uj)dt + V2Kdw{t), x{0) = 0, 

where x{t) is the position of the particle at time t, k > the molecular diffusivity, w{t) 
the standard Brownian motion and V(x, t, a;) a time-stationary, space-homogeneous, incom- 
pressible velocity field. Here u denotes an element of an ensemble of random flows. 

We are concerned with the long time, large scale behavior of the displacement x{t). To 
this end, we study the scaling hmit 

x%t) =ex{t/6'^''), ase^O, (1) 

with suitable g > 0. The scaling exponent q characterizes the time scale associated with 
transport on the spatial observation scale 1/e. The equation for the rescaled displacement 
(0) becomes 

dx%t) = e^-^''y{x'{t)/e, t/e^'')dt + e'^-''V2^dw{t). (2) 

When molecular diffusion is evidently negligible, we set k = to simplify the equation 

dx'{t) = e'-^''V{x%t)/e, t/e^'')dt. (3) 

The effect of molecular diffusion is discussed where the issue arises and in the concluding 
remarks. 

Motivated by existing diffusion limit theorems for steady flows with finite-range, spatial 
correlations ([^,13), on one hand, and those for temporally mixing flows with long-range, 
spatial correlations (0, |]T5|, |31|), on the other hand, we consider turbulent transport in 



a class of random flows with power-law spectra parametrized by a,/5,7 (see Section 2 for 
details). 

Roughly speaking, the velocity field V is time-stationary, space-homogeneous and Gaus- 
sian. Its two-point correlation functions R = [i?jj], i?jj(x, t) = (?7j(-, ■)Uj{- -|- x, ■ -|- t)), are 
given by the Fourier transform Rij{'yc, t) = J e*'^'^-Rij(k, t)d]<. with 

R.,{k,t) = p{\k\^^t)S{k) {6,j - hk,\k\-^) \k\^-'' (4) 

where p is the time correlation (relaxation) function and £ the energy (shell) spectrum given 
by a power-law 

£(k) = ^o|k|^"^", ^o>0. 

Here (■) denotes the ensemble average. The factor {5ij — kikj\\i.\^'^) ensures that the flow 
is divergence-free. If p is an exponential function, e"'^"'''' *, then the velocity field is an 
Ornstein-Uhlenbeck process which is Markovian. 

Note that the spectrum is not integrable near |k| = oo or |k| = for a < 1 or a > 1, 
respectively. The infrared divergence (small |k|) of the integral of velocity energy spectrum 
indicates non-homogeneous velocity and thus violates the space homogeneity assumption, 
whereas the ultraviolet divergence (large |k|) of the integral would make the velocity a 



generalized, rather than ordinary, function (i.e. a distribution). To remove divergence in the 
spectrum, we introduce an ultraviolet cutoff 



and an infrared cutoff 



£{k) = 0, |k| > K, for a < 1 (5) 



f (k) = 0, |k| < (5 < 1, for a > 1. (6) 



In the case of a < 1, the energy containing scale is at the ultraviolet cutoff; in the case of 
a > 1, the energy containing scale is at the infrared cutoff. It is convenient to write the 
cut-off energy spectrum as 

£(k)=i?o|k|^-2"J(|k|) (7) 

where /(|k|) is the characteristic function of [0,/^], for « < 1, of [5, oo) for a > 1, and of 
[(5, i^], for a = 1 (see |^). When we study the effect of an infrared cutoff, we will take 
an infrared cutoff 6 = e'^ > 0. Ultraviolet cutoffs play less prominent role than do infrared 
cutoffs in the scaling limit. 

If the infrared cutoff 5 > is fixed, independent of e, then the flow is mixing in time (i.e. 
correlation time is uniformly bounded, independent of wave number) and, consequently, the 
scaling in (|1]) is diffusive, q = 1, and the limit is a Brownian motion ([0). In the case of the 
Kolmogorov-Obukhov spectrum (a = 4/3, /5 = 1/3, see [^), [S,K] represents the inertial 
range where K~^ is the dissipation length and 6^^ is the integral length. In general, the 
infrared cutoff is determined by the scale of external forcing and the size of physical domain. 
By letting 6 change with e, as 6 = e'^,7 > 0, we vary the spatial scale of observation 1/e in 
relation to, e.g., the size of physical domain. 

If the scaling limit exists, statistically independent of the initial point, and has station- 
ary increments, then the transport process is said to be homogenized, and the up-scaling, 
or coarse-graining, procedure represented by (|l]) is called homogenization. The scaling is 
diffusive if g = 1, super diffusive if g < 1, subdiffusive if g > 1. Sub- and super-diffusions are 
called anomalous diffusion. 

The limit Z(t) may be Gaussian or non-Gaussian, Markovian or non-Markovian, even 
if the velocity field is Gaussian and Markovian. In general Z(t) has stationary increments 
as does x^(t) ([^). If Z(t) is self-similar and Gaussian then it can be characterized by a 
unique Hurst exponent H in its autocovariance function 

Cov(Z(ti), Z(t2)) = ic {|ti|2^ + |t2r^ - |ti - t2r} , 0<H<1 (8) 

where C is the variance of Z(l). H = 1 corresponds to a regular (or smooth) motion; H = 1/2 
a Brownian motion, B(t). Any other H corresponds to a fractional Brownian motion (FBM), 
Bi/(t), which, after normalization, can be represented as 



B 



H 



(t) = r {\t - t'l^-i/2 _ |t'|^-V2)rfB(t') + /* |t - t'|^-^/2^B(t'), < if < 1, (9) 

J-oo JO 



as introduced in [^. Eq. (^ defines the only mean-zero, mean-square, continuous, Gaus- 



sian process that is self-similar (or self-affine), with the Hurst exponent if, has stationary 
increments, and satisfies Bj:/(0) = (see ^l|). FBMs found in the present study are all 



persistent in the sense that H > 1/2. It is worth noting that, for a > 1 with the critical 
cutoff (7 = 7c), the hmit process is not self-similar (cf. Regimes II', III' and IV'). 



Non-Markovian limits are related to non-local homogenization (Pj, |[TT|], [^, |^4[). Pre- 
viously nonlocal homogenization has been shown to arise as a result of fast oscillation, rather 
than of a scaling limit. 

If Z(t) is not Gaussian, then there may be a hierarchy of Hurst exponents corresponding 
to higher moments of the process. When the sequence of Hurst exponents diverges as the 
order of moment increases, the limit is intermittent. Intermittency effect may also manifest 
as multiple scaling exponents. We do not consider the problem of intermittency here. 

In this paper, we do not address directly the question of existence and uniqueness of the 
scaling limit. Rather, we assume the existence and uniqueness of a nontrivial scaling limit, 
and seek to identify the scaling exponent and the second order Hurst exponent (the Hurst 
exponent, in case of a Gaussian limit). In doing so, we point out relevant mathematical 
results that exist, or can be or are yet to be proved. We try to present a simple, coherent 
physical picture of the whole phase diagrams. To enhance our case, we often analyze the 
problems from several different perspectives. 

The exponent q characterizes the time scale associated with transport observed on the 
space scale 1/e; the exponent H characterizes the time correlation property of successive 
increments on the observation scale (and, therefore, the roughness of the limiting sample 
paths). Naturally we ask if the dimensionally correct relation 

H = l/(2g) (10) 

holds? When (|T0|) holds the limit process is invariant under the same scaling transformation 
(0). It turns out that relation ([T0| ) generally holds for a < 1 but fails for a > 1. If an 
additional infrared cutoff is made in the case of a < 1 and if the cutoff is removed faster than 
some critical wave number, kc, then ([T0|) does not hold. In these situations, the inequality 

H < l/(2g) (11) 

is in place of (p!0|). The inequality (jTl|) is due to the fact that 2H characterizes the covari- 
ances, whereas 1/q characterizes the variances, of successive increments of turbulent motion 
on the observation scale. 

In general the exponents q, H depend on the parameters a, f3 and the cutoff 6 = e'' and 
can be expressed exphcitly as functions of a, /5, 7. Here it may be helpful to draw analogy 
to critical phenomena in statistical physics: we think of a,/3, 7 are order parameters and 
the scaling limit e ^ as thermodynamic limit and the exponents q, H given by formulas 
of a, P, 7 as phases. The phase diagram divides the space of order parameters a, P, 7 into 
regions associated with different formulas for q, H. Our results are summarized in Figures 
1, 2 , 3 and 4. Since there are three parameters, the full phase diagram is three dimensional. 
To simplify the presentation, we choose to portray the full diagram as two two-dimensional 
diagrams, one for supercritical and the other for subcritical infrared cutoffs. 

Note also that the phase diagrams are different from those in statistical mechanics in that 
our phases are continuum, not discrete: except for the diffusive regime of a + /3 < 1 where 
H = 1/2, q = 1, H,q change from point to point, continuously or discontinuously. But their 
functional forms in relation to a,/9, 7 are discrete and divided by phase boundaries. 



Phase diagram was first used by Avellaneda and Majda [0, to present scaling limits 
of turbulent transport in anisotropic, stratified flows of the form V(x, t) = {v{x2, t),0), with 
X = (xi,a;2). A different diagram for the same shear- layer flows was rigorously obtained by 



Zhang and Glimm [^ using a different approach. In the current paper we consider isotropic 
turbulent flows and the results are different from those for anisotropic flows. Also, we do 
not attempt to derive the results rigorously here. Often we refer to existing theorems to 
indicate how in principle results may be proved, subject to technical modification, and to 
support the physical arguments invoked; they are not intended to be mathematical proofs. 
The proofs of many of the results are very technical and will be published elsewhere. 

The effect of an infrared cutoff depends on whether the cutoff is subcritical or supercritical. 
For a < 1, a supercritical cutoff, 7 > 7c = max{l, l/(a + 2/3 — 1)}, does not affect the 
scaling law. For a > 1, because the infrared cutoff corresponds to the energy-containing 
scale, the scaling limit is dominated by the infrared cutoff. 

The supercritical diagram includes: 

• Regime I:a-|-/3<lorQ;<0. The scaling is diffusive, g = 1, and the limit is a 
Brownian motion, H = 1/2. 

• Regime II: a + /3 > 1, a + 2/3 < 2, a < 1, 7 > l/(a + 2/3-1). A FBM regime with the 
space-freezing property that velocity dependence on space is negligible. The scaling 
is superdiffusive, q = (3 /{a + 2/3 — 1), and the limit is a fractional Brownian motion, 
H=l/{2q). 

• Regime III: a + 2/3 > 2, < a; < 1,7 > 1. A FBM regime with the time-freezing 
property that velocity dependence on time is negligible. The scaling is superdiffusive, 
g = 1 — a/2, and the limit is a fractional Brownian motion, H = l/(2g). 

• Regime IV: 1 < a < 2, 7 > max{l, l/(a + 2/3 — 1)}. A regular (or smooth) motion 
regime with both the space-freezing and the time-freezing properties. The scaling is 
superdiffusive, g = (1 + 7)/2 — 70/2, and the limit is regular {H = 1). 

The relation (^) is satisfied in all but Regime IV. 

In the case of subcritical cutoffs, 7 < 7c = max{l, (a + 2/3 — 1)~^}, the number of 
regimes shrinks as the significance of low wave numbers is reduced: part of Regime IV merges 
with Regime II, and part of Regime IV merges with Regime III. The scaling exponent now 
depends on the cutoff exponent 7 explicitly. The limit is universally a Brownian motion 
across all regimes. 

The subcritical regime includes: 

• Regime I remains intact. 

• Regime V: a + (3 > 1, a + 2/3 < 2, 7 < l/(a + 2/3 — 1). Velocity decorrelation in time 
dominates the transport. The scaling is superdiffusive, g = 1 + 7 — 7(0 + /?). 

• Regime VI: a + 2/3 > 2, a > 0,7 < 1. Velocity decorrelation in space dominates the 
transport. The scaling is superdiffusive, g = 1 — 70/2. 

Finally, there are three regimes associated with critical cutoffs for which the limit process 
is not self-similar and, thus, the Hurst exponent is not well-defined. 
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• Regime IP: a + p > l,a + 2p < 2,0 < a < 1 with 7 = (a + 2/? - 1)"^ 

• Regime III': a + 2(3 > 2,0 < a < 1 with 7 = 1. 

• Regime IV': 1 < a < 1 + I/7 with 7 = max {l,{a + (3 - 1)"^}. 

Part of Regimes V and VI was first studied by Avellaneda and Majda |^ (see also 
!l||). The main difference in assumption and setup between this work and is that they 
considered a partial diagram (0 < /? < 1/2, < a < 2) with an infrared cutoff 7 = 1 < 7c 



(see also |^). Figure 4 is a generalization of theirs. The phase diagram of ||23[ was obtained 
entirely by certain scaling arguments, and is restricted to two dimensions. 

By contrast, our main findings are: (i) the transport effect of the sampling drift and 
related critical wave number, which are introduced for the first time, (ii) fractional Brownian 
motion limit as a result of the critical wave numbers, (iii) the effect of infrared cutoffs, (iv) 
the formulation of cutoff dependent eddy diffusivity and its associated variational principle 
without molecular diffusion. The variational principle gives a useful bound for the eddy 
diffusivity. 

The organization of the paper is as follows. In Section 2, we define the three-parameter 
family of Gaussian flows, whose transport properties are discussed in subsequent sections. 
In Section 3, we introduce the notions of sampling drift, critical wave numbers and eddy 
diffusivity. We also formulate variational principles that lead to general bounds for the 
cutoff dependent eddy diffusivity in terms of a fractional vector potential of the velocity 
field. Since the transition from ultraviolet to infrared cutoff in velocity occurs at a = 1, we 
divide the discussion accordingly into two cases: a < 1 and a > 1. We consider the case 
a < 1 in Section 4 and the case a > 1 in Section 5. We conclude with various remarks 
in Section 6. In Appendix we derive a variational principle for the cutoff dependent eddy 
diffusivity, without the presence of molecular diffusion. 

2 Random velocity field 

In this section, we describe some mathematical properties of the random velocity fields 
considered in this paper. 

The most important property is stationarity in time and homogeneity in space (space-time 
stationarity for short), without which homogenization is unlikely to hold. It should be noted 
that, when formulated in a general, abstract framework as we will do momentarily, space-time 
stationarity encompasses space-time periodicity, quasi-periodicity and almost periodicity as 
well as random stationarity. This abstract formulation is also handy for formulating the 
variational principle for the eddy diffusivity (Section 3.2). Elsewhere, the paper can be 
understood without referring to the abstract formulation. 

The variational principle in the absence of molecular diffusivity also uses explicitly the 
Markov property of the flow and the associated generator. A key turbulent diffusion theorem 
cited in the discussion of the diffusive regime (Section 4.1) was proved for certain Markovian 
velocity fields. For Markovian flows, the mixing property conveniently corresponds to the 
spectral gap of the generator. Elsewhere, the Markov property is not used explicitly and 
probably not needed. 



Since we only use the (energy) spectral density explicitly in presenting the phase dia- 
grams, it is safe to assume that the velocity fields are Gaussian. In particular, the Gaussian 
property is essential in the fractional-Brownian- motion regimes (II and III). Elsewhere, the 
Gaussian property is probably not important. 

Let us begin with the abstract formulation of space-time stationarity, upon which we 
will define the Gaussian and Markov properties. Let ft be the space of steady, space- 
homogeneous velocity fields and let P be a probability measure on ft. Homogeneity in space 
can be canonically described by the invariance of the distribution P under the group of 
translations {Tx}xeR<* acting on ft. We further assume that P is ergodic with respect to 
{Tx}xe/?<* ill the sense that the only invariant, measurable functions on ft under {Tx}xe/?d 
are constants. The measure P dictates the correlation of the velocity field in space, and, in 
case of Gaussian velocity fields, is determined by the energy spectrum. 

Alternatively, we think of fi as the ensemble of elements u, representing the randomness 
of the velocity field, which is distributed according to the measure P. A (prototypical) 
random velocity field is a vector- valued, random variable (i.e., a function on f2), denoted by 
V(co'). The realization or the sample of the (time independent) velocity field, V(x, cu), is the 
translate of V(ci;) on O, i.e., V(x, a;) = V(txI-i^). Since the measure P is invariant under 
the translations, the resulting velocity fields are space-homogeneous. We assume that V has 
zero mean 

(V) = 

and zero divergence 

V-V = 0, V = idi,d2,...,dd). 

Partial derivative di is the infinitesimal generator of the subgroup of translation {r^^.j^^.g/j. 
The Laplacian A := V ■ V is defined as usual. As before, (■) denotes the ensemble average 
with respect to the distribution P. 

The time dependence of the velocity field is then introduced as a time-stationary stochas- 
tic process, uj{t), on the space O, which preserves the measure P. In other words, P is an 
invariant measure of the process uj(t). The realization of time dependent velocity field is 
then given by 

V(x, t, uj) = Y{T^uj{t)), uj{0) = u. 

In this formulation, the temporal properties are conveniently separated from the spatial 
properties of the velocity field. Additional structures such as Gaussianity and Markovianity 
can be added on by imposing corresponding properties on P and uj{t). The space O is 
usually infinite dimensional in suitable coordinates such as Fourier modes. This formulation 
is sufficiently general to describe periodic, quasi-periodic, almost periodic as well as random 
homogeneous velocity fields (see, e.g., [p!0|l ). 

We think of a Markovian velocity field as a sample path in f2 of a Markov process uj{t). 
A Markovian, Gaussian velocity field corresponds to an exponential time correlation function 
p in (H) and admits the spectral representation 



V(x,t)= / e'2"''-"V(dk, t) 



_Rd 



where the stochastic measure V((ik, t) is an Ornstein-Uhlenbeck process 
dtV{dk,t) = -ao\k\^^V{dk,t)dt + \kfS^/\k) (l-k^k\k\-^Y^'^ \k\'^^-'^^/^W{dk,dt) (12) 



and can be conveniently expressed in terms of a Gaussian white noise W((ik, ds) 
t 



~oo 



The Ornstein-Uhlenbeck process ([T^) has an invariant measure P that is a Gaussian 
distribution with zero mean and the variance matrix R = [Rij] given by (1). Then the 
exponential relaxation function corresponds to a generator A of the form 

^=(-^AfA, P>0 (13) 

where ^o is the generator of the process 

dtVo{dk,t) = -aoVo{dk,t)dt + S^/\k) (l - k ® kjkj^^)^^^ |k|(i-'^)/2w(t/k, rft). (14) 

The operator Aq is symmetric with respect to the measure P and commutes with the trans- 
lation Tx, Vx G W^. As the process (|1^) is a time change of (p!2D and different wave numbers 



are independent, the measure P remains invariant with respect to (|T^. Also, because the 
time correlation function for (|T^) is exponential with an exponent Oq uniformly bounded 
above zero, ^o has a spectral gap 

- {Aoff) > aoiP), ao > 0, (15) 

for all functions /, (/) = 0, in the domain of ^o- 

The motion in this temporally stationary, Markovian flow is also a temporally stationary, 
Markov process whose generator is 

C = A + V-V (16) 

when molecular diffusion is absent, and is 

C = A + kA + V-V (17) 



when molecular diffusion is present {^l 

Now we make an observation which will be used later in assessing the role of molecular 
diffusion. The generator ([T7|) in conjunction with ([T3|) and ([T5| ) suggests that the presense 
of molecular diffusion introduces a mechanism of generating a Lagragian correlation in time 
comparable to /3 = 1 in the Eulerian correlation in time. For (3 < 1, the generator A 
dominates over kA for low wave numbers and, if a fixed ultraviolet cutoff is also present, are 
also comparable to kA for the other wave numbers. Thus, the effect of molecular diffusion 
is negligible for (3 < 1 and a < 1 in the limit of high Peclet number (k ^ 0). 

In the sequel we shall use the notation of the fractional gradiant of order /3, 

V/3 := (_A)(^-l)/2y. 



3 Transport properties of various wave numbers 

To study motion in a flow with a power-law energy spectrum over a wide range of scales, it is 
convenient to decompose the energy spectrum into the the sampling drift and the fluctuating 
velocity field, and to consider separately their distinctive transport properties. The relation 
between the sampling drift and the fluctuating velocity field is like that between a mean flow 
and the fluctuation. 

3.1 Sampling drift and critical wave numbers 

For each realization of random velocity field there is a nonzero sampling drift due to random 
fluctuation, depending on the scale of observation. 

The volume-averaged flow on the observation scale consists of spatially non-fluctuating 
wave numbers on the observation scale, namely, all |k| ~ £. The the volume-averaged flow 
comprises three kinds of wave numbers: supercritical, critical and subcritical wave numbers 
depending on their variations in time on the observation scale. Critical and supercritical 
wave numbers compose the sampling drift. 

The supercritical wave numbers are effectively steady in the sense that their correlation 
times are much larger than the time scale of observation, |k|^^^ ^ e^'^'^, thus, satisfy 

|k| <min{e''/'^,e}. (18) 

As these wave numbers are temporally as well as spatially uniform, they behave like a 
constant drift on the observation scale and transport particles ballistically. Among them, 
we pay special attention to those wave numbers that, on their own correlation time scales, 
transport particles over a distance larger than the observation scale 

|k|i-"|k|-2^ > 1/e, (19) 

since wave numbers of order |k| have an amplitude of the order 



/ci|k|<|k'|<C2|k| 



\ 1/2 

£{k')d\k'\ ~ |k|i-", |k| < 1. (20) 



Note that, for (|19]) to define a non-empty set of low wave numbers, we need a + 2f3 > 1. 



For a + 2P < 1, the supercritical wave numbers do not contribute to the transport on the 
observation scale and are negligible asymptotically. In any case, the transport effect of the 
insignificant supercritical wave numbers are by nature negligible. 

Since we do not know the scaling exponent q ahead of time, we define the critical wave 
numbers to be the boundary of those significant supercritical wave numbers. Thus, the 
critical wave numbers are of the order kc = e'^", with 

/W ^oR ^\-n j {a + 2(3-l)-\ ioTl<a + 2(3<2 . 

7. = max {1, (a + 2/3-1) } = | ^^ for a + 2/3 > 2. ^^1) 

By (p]), for a + 2/3 > 2, the the sampling drift is identical to the volume-averaged flow. 
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Insignificant supercritical wave numbers occur when the following conditions are satisfied: 
7c > q/ (3 > 1. This leads immediately to a + 2/? < 2 and /3 < 1. The latter follows from 
g < 1 (see Section 4.1). As we will see later, this can only happen, in part of Regime I (with 
q=l) defined by a + /3 < 1, a + 2/? < 2, /? < 1. 

The critical wave numbers have long-range correlation in time or in space and dominate 
the transport in the fractional-Brownian-motion regimes (Regimes II and III). The suhcritical 
wave numbers are either temporally fluctuating, |k| ^ e*^/^, or spatially fluctuating |k| ^ e. 
Effectively, the subcritical wave numbers can be defined by |k| ^ k,, and, by definition, 
include the insignificant supercritical wave numbers. 

Denote by q the sampling drift on scale 1/e. It has an amplitude of the order 

([ crtwinV^' / I A;!-" -51-" I, fora^l ,_. 

iW/^'^^'"'^"; -||logfc,-log5|, fora = l. (22) 

Since the critical wave numbers dominate the sampling drift for a < 1, c^ has a long-range 
correlation in space or time on the observation scale, so its transport effect is not ballistic. 
For a > 1, Cj is effectively frozen in time and its transport effect is ballistic. 

Infrared cutoffs are classified accordingly: 5 = e'^ is critical if 7 = 7c, supercritical if 
7 < 7c, and suhcritical if 7 > 7c. We call 7c the critical exponent. 

From (|18]) and ([T9|) , we have the simple inequality for the scaling exponent 



/ P/{a + 2/5-1), for 1< a + 2/3 < 2 
^ - \ /?, for a + 2/3 > 2 ^^"^^ 

The equality in (p3D is admissible because asymptotics is a continuum and can not be fully 
resolved by power-laws. Here we restrict our attention to the power-law part of scaling 
behaviors. 

For transport effect, besides the line a + 2/3 = 2, the line a = 1 is also important for 
the following reasons. For a < 1, the sampling drift is dominated by the critical wave 
numbers, whereas, for a > 1, the sampling drift is dominated by wave numbers nearby the 
infrared cutoff. Moreover, in the case of a > 1, the infrared cutoff corresponds to the energy 
containing scale and, therefore, dominates the the transport as well as the flow. As a result, 
scaling laws of transport for a > 1 are in general (infrared) cutoff dependent. The limit 
processes in the case of a supercritical cutoff, however, are always regular motions [H = 1, 
Regime IV) as the effectively constant drift dominates the transport. 

Based on the supercritical wave numbers alone, the exit time r (out of a ball of radius 
1/e) for a < 1 can be estimated by 

, _ r .-W(«+2,-i)^ for 1 < a + 2/? < 2 
^^"^ /^ ~ 1 £-2+", for a + 2p>2. ^^^^ 



It is easy to see that, the (asymptotic) equality in ( p4D is achieved when the combined effect 
of the supercritical and the critical wave numbers is considered since, for a < 1, the critical 
wave numbers are much stronger than the supercritical wave numbers in magnitude. For 
a > 1, however, the transport is dominated by the wave numbers |k| ~ 5. So we have 

r < S^'-^/e = £-i-T+"T (with 6 = e^) (25) 
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in the case of a > 1 . 

As we will show by the variational method in Section 3.3 that the critical wave numbers 
dominate the transport in Regimes II and III. In the case of a < 1, the exponent for 
a + 2(3 < 2 is less than or equal to 2 (i.e., 2/?/(a + 2/? - 1) < 2) only ii a + (3 > 1; for 
a + 2(3 > 2, the exponent is less than or equal to 2 (i.e., 2 — a < 2) only if a > 0. The former 
defines Regime II; the latter defines Regime III. In the case of a > 1, any nonnegative 7 
leads to 1 + 7 — 07 < 2 (the scaling is superballistic for 7 > 1). In the remaining region 
(Regime I: a + /? < 1 or a < 0), the supercritical wave numbers are negligible since the 
transport effect of the fluctuating wave numbers is at least diffusive as we will see later. 
Equating the exponent with 2g, we have, from (p4|), the scaling exponents for Regimes II, 
III (see Section 4), and, from eq. (pSf), the scaling exponent for Regime IV (see Section 5), 
both with supercritical cutoffs, 7 > 7c. 

In the regimes where the critical wave numbers have a leading effect, the scaling limit is 
a fractional Brownian motion (Regimes II and III). Fractional Brownian motions arise as a 
result of long-range correlation of the critical wave numbers. 

If the infrared cutoff is subcritical, i.e., 5 ^ fcc, wave numbers of the spectrum are 
either temporally or spatially fluctuating. Contrary to the fractional Brownian motion limit 
caused by the critical sampling drift, the limit is always a Brownian motion. But the scaling 
exponent may be superdiffusive due to low wave numbers in the vicinity of the cutoff. 

3.2 Subcritical wave numbers: eddy difFusivity 

To study the effect of subcritical, or fluctuating, wave numbers on transport, we think of 
turbulent motion as a superposition of a mean flow (i.e. c^), and the fluctuating flow, 
following a spectral discretization. 

We propose that the fluctuating wave numbers give rise to a fluctuating motion, on top 
of the mean flow, on the observation scale, and this fluctuating motion can be characterized 
by a notion of scale dependent eddy diffusivity introduced below. We then formulate two 
variational principles and use them to obtain general upper bounds for the (scale-dependent) 
eddy diffusivity. 

Spectral discretization is motivated by a standard result of the ergodic theory for sta- 
tionary processes that a stationary process is a limit of periodic processes (see [|10|)- We will 



use the periodic approximation in two different ways: In the first, we consider the periodic 
approximation in the space variables only and work with a subspace of (ri,P), the space 
(fi*-"^, p(")) of time-independent, space-periodic velocity fields with period cell [0,n]°' (see 
discussion below). In this approach, time randomness in the velocity field is represented as a 
Markov process on (f2''"\ P^"')). In the second approach, we work with a sequence of space- 
time periodic fields with the (normalized) Lebesgue measure as the probability distribution 
on the space-time period cells as stated in the following lemma. 

Lemma 1 Let u be a stationary process. Then there exists a sequence of periodic processes 
ujn of period in ^ 00 in each variable, such that, the probability measure Pn obtained as the 
distribution of r^^uon where x is random and distributed uniformly on the period cell [0,£„]'^ 
converges weakly to the distribution of uo as n ^ 00. 



12 



(See, for instance, [37| for a proof). We emphasize that spectral discretization is only a con- 
venience for the formulation of the variational principles; it is neither essential nor necessary. 

We now formulate the first approach more specifically. A spectrally discretized flow can be 
written as a sum Ce + V^^'"^ where c^ is the sampling drift (see Section 3.1 and eq. ([2^)), and 
yie,n) jg ^]^g spatially periodic version of the fluctuating velocity field with a discrete spectrum 
k G Z'^/n,max{kc,6} < |k| < ii^ and the amplitude (I-k(g)k/|kp) J/|k|<|k/|<|k|+i/n^(k')f^|k'|. 
The mesh size 1/n should tend to zero sufficiently fast, as e — >■ 0, to approximate the trans- 
port effect of the original fluctuating flow in view of the above lemma. 

Equivalently, we replace the spectral measure V((ip, t) by the discrete measure V*^'^'"^(k, t)(5k,p, 
Vk G Z'^/n,max{kc,5} <^ \k\ < K with V(^'")(k,t) satisfying 



ciiV(^'")(k,t) = -ao|kp^V(^'")(k,t)dt+|k|^J / 



k|<|k'|<|k|+l/n 



S(k')d\k'\(l-k^k\k\ 



oxl/2 



dtW{k,t) 



(26) 

where W(k, t),Vk G Z'^/n,Taeix{kc,S} <^ \k\ < K are independent standard Brownian 
motions. As discussed in Section 3.1, the sampling drift c^ is steady for a + 2/5 > 2 or a > 1; 
it is unsteady for a -|- 2/? < 2, a < 1. 

The time-stationary, space-periodic field V^^'"'^(x, t, c<j„), cu G ri'"-* is a Markovian fiow and 
can be represented as a translate, V('^'")(x, t, cun) = V'^^'")(x, ujn{t)), of steady, space-periodic 
field V^^'"^(x, u;„) where a;„(t),'u;„(0) = is a Markov process on Ct^^'. As usual, we write 
ujn explicitly only to emphasize its role. 

For fixed e,n, the displacement, x{t), in the periodic fiow, c^ + V'^^'"\ is the sum of a 
mean motion, jQC^{s)ds, and the fluctuation, x(t) — jQC^{s)ds. After a proper rescaling 
t — > A^t, X -^ Ax, A -^ CX3, the fluctuation converges to a Brownian motion by a turbulent 
diffusion theorem for mixing flows ([|T^). Let ^(^'") be the generator for c^(t) + V'^^''^)(x, t). 

{e,n) 



The diffusion coefficients, D^^ 

random, space-periodic solution Xj " 

of the abstract cell problem (cf. (0), see also [p!5| ) 



of the limiting Brownian motion are determined from the 
i.e. Xj can be viewed as a function defined on fi^"-*) 






V/''"\ in rj("), V^,j (27) 



D. 



{e,n) 



^ i ((\/f '")xf '"^)n + (V^/^'"^xf"^«) 

= (V^xf '"^ ■ 4^'"V^X?'"V, V^,j. (28) 

with the periodic boundary condition, where (■)„ is the average with respect to P^'^\ Here 
we have used the following identity 



(e,n) 



([(Ce+V(^'"))-Vx 

= v-((Q + v(^''^))xf'")xf'"^) 



x\ 



(^'"))„+([(c,+v(^'"))-vxf"^; 





Xi ) 



(29) 



which follows from the incompressibility of c^ + V^^'"-* and the space-homogeneity of ((c^ + 



Y(£,n)^^ls-"J^te>"i^^_ 
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The problem (^) is well-posed and has a unique solution, up to a constant (which does 
not affect (pH])). In Appendix, we derive the minimum principle 

D(^'")(e) := D^^-'^^e "6 = 1^ {-{A^''^^ff)n - (^^''"VT)n} (30) 

with the space-periodic functions /', / related by 

^fe")j' + (c, + V("'"))-V/ + V("'")-e = 0, in O^"). (31) 

It should be noted that the explicit form of the generator is not used for the variational 
formulation. 

In the limit n ^ oo, the abstract cell problem (P7])-(P^) becomes 



Cxf + V,^"^ = 0, infl. (32) 



We also have 



Df.> = lim D„ 



^'^ - limDf^'") 

5 {lyr^') + (vf 'xf )) 

(Vxf'-AVxf), Vi,j, (33) 



following from (|2^). It is clear from ( P^D that the matrix D*-^^ = [-D^j ] is symmetric and 
positive-definite. We think of D*^"^) as a measure of turbulent dispersion caused by eddies 
composed of subcritical wave numbers in interaction with the sampling drift. We call it the 
eddy dijfusivity. If the increments of the fluctuation of particle motion have divergent step 
sizes as e — > 0, then the eddy diffusivity is cutoff-dependent. Eq. ( |5BD indicates the right 
solution space for (|3^): L|(f2), the space of functions with homogeneous, square integrable 
fractional gradient of order (3. 

From another perspective, the variance of the fluctuation x{t) — jQC^{s)ds after the 
rescaling t -^ A^t, x —>■ Ax, A ^ oo can be expressed as the time integral 



2j^D',^\s)ds 



'0 

of the Lagrangian velocity autocorrelation 

Dt-'is) = \ £^' {{Vl^\x{s), s)vl^\x{s'),s') + {Vj'\x{s), s)vl^\x{s') , s'))) ds' . (34) 

Because the Lagrangian velocity V*^^)(a;(t), t) is a stationary Markov process (|^,|T^), ( [5^ ) 
can be rewritten as 

Df{s) = I £'' {{Vl'\0,0)v}'\xis'),s')) + {V}'\0,0)vl'\xis%s'))} ds' 
= I jf^ {{Vl^^ exp {Cs')vl'^) + {Vl'^ exp iCs')vl'^)} ds'. 
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In the limit A -^ oo, Df^ (s) tends to the following expressions 



(v^'x^') + ^^^) 



D 



(e) 



(35) 



where xf is the solution of (|32D . 

When molecular diffusion is present, we denote the eddy diffusivity by D^'^'"'^ As before, 
j^{e,n) ^g^^ i^g characterized variationally by adding the terms, k, fi;(V/ ■ V/)n, f^i^ f ■ ^ f')n 
to ( PU[ ) and a Laplacian term, kA/', to (PT|). 

We turn to the second approach of space-time periodic approximation. Let \'^'^'^'^\x,t) 
be the approximating sequence of space-time periodic fields, as stated in Lemma 0, with 
increasing space period n and time period A, for the velocity field V'-^-* with a subcritical 
cutoff. We work with the space-time period cell problem in which time randomness in the 
Lagrangian dynamics is absent. To formulate a variational principle in this case, we need to 
reinstate the molecular diffusion here. 

For fixed A, n, the effective diffusivity, D^f'"''^-', in the flow V''^'"''^-*(a;, t) exists and can be 
given as 

D^"'^)(e) =mi^ f'^ f Kil + V/ • V/ + Vf ■ V/') rfx dt (36) 

/ A JO W^ J[0,n]'* 

where /, /' are both temporally and spatially periodic with the period A, n, respectively, and 
are related to each other by the following equation 

dl 

dt 

1% and ||19|). The eddy diffusivity D^*^) in the presence of molecular diffusion is the 
large scale limit of D\^'^'^\ i.e., 



/.A/' 



[Ce + V("'"'^)) ■ V/ - V 



(e,n,A) 



see 



D(^) 



lim D^f'"'^). 

n,A— ^cxD 



(37) 



The variational principle (p^)-(p7|) is more useful than (|30|)-(|3T|) when the Laplacian 
in the generator ([T7| ) dominates A for low wave numbers as in the case of /? > 1. Another 
advantage for working with the space-time periodic setting is that a dual variational principle 
can be formulated for the inverse of D^^'"^'^^ and can be used to obtain the lower bound for 



VI 



3.3 Variational bounds: fractional vector potential 

3.3.1 Case 1: Supercritical cutoff 

When the sampling drift is present, i.e., 7 > 7c, we show by the variational principles the 
following upper bounds on the growth rate of the eddy diffusivity 



D 



(e) 



< C, ior a + /3 <1 or a <0 



Dff « 



^27e(l-a-/3) 
log(l/£>),' 



for a + /3 > 1 
ioY a + P = 1 
for a > 



(38) 



log(l/£>), for a = 
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Vi,j, for some constant C > 0. Note that e " is a better bound than e'^"'"^^ " '^^ for 
a + 2(3>2. 

Take the trivial trial function / = in (^) and eliminate the first term in the functional. 
We calculate the second term in ( pOD by studying the equation 

^fe")/' + V("'")-e = (39) 

(cf. (0)). Consider the fractional vector potential (the fractional stream function) in three 

^/3 



dimensions (in two dimensions) HL^'" of order /3 defined by 



via the Fourier transform. This means that HL^ = lim„^oo HL^'" has the energy spectrum 

with a subcritical infrared cutoff. The usual vector potential and the stream function cor- 
respond to /? = 1. What is significant is that, for a + /3 < 1, ( |^ is integrable near k = 
uniformly as the the infrared cutoff is removed, and thus defines a homogeneous, square inte- 
grable H/3 that is cutoff- independent. For a + /3 > 1, ( |iT| ) is not square integrable uniformly 
as the the infrared cutoff is removed, and the second moment of H'g grows like 

iinfl^) < £27.(i-a-/3)^ for « + /5 > 1 (42) 

and (|hJ^|2) < log(l/e^O, for a + /5 = 1 as e -^ 0. 
In terms of Up"'' , (^Uj) becomes 

^ie,n)j, ^ (_a)^/2h^^'") . e = 0. (43) 

A straightforward energy estimate for (^) gives 

- (^(^'")/7')n = ((-A)^/2Hj'").e/% (44) 



r(e,n) 



< V(|Hr ■eP)W(l(-^)"/'/'l')n (45) 



< y(|<")-eP)„^l(-^(^.")/'f)n. (46) 



Therefore 



-(^(^'")/r)„ <-(!<") ■ep)„ (47) 

ao 

which, in the limit n -^ cx), is much less than 527c(i-a-/3) f^j. g^all e. 

It is worth noting that the right side of (|^) is, up to a factor independent of e, what 
one gets in replacing the Lagrangian autocorrelation in ( p4D by the Eulerian autocorrelation 

I^^R^\0,s)ds (48) 
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which is called the (Eulerian) Taylor-Kubo formula, used extensively in the literature to 
approximate the eddy diffusivity since Taylor's work ([0). The physical significance of the 
bound (0) is that the eddy diffusivity Df^' of the fluctuation is bounded, as the infrared cutoff 
is removed, by a constant times the Eulerian Taylor-Kubo formula ^4^); the eddy diffusivity 
may be much smaller than (^) due to the spatial decorrelation of velocity. 

A different upper bound for the eddy diffusivity the can be obtained by using the second 
variational principle (|36|). First we note that the eddy diffusivity in the presence of molecular 
diffusion would be enhanced if we freeze the time variable of the velocity field V*-^-*. This 
can be easily seen as follows. Let D^^'"^ be the eddy diffusivity for the frozen velocity field 
V(^''^)(x,0) 

D^:'^\e) := D(f'")e ■ e = ini -J n{l + Vf-Vf + Vf ■ Vf) dx (49) 

where /, /' are spatially periodic with period cell [0,n]'^ and are related by 

nAf = -(c, + V("'")(a;, 0)) • V/ - V("'")(a;, 0) ■ e. 

Since time independent trial functions / are admissible in (|36D , (^91) is larger than D^'^'"'^\ VA > 
0, given by (|5B|). 

Using the trivial trial function / = in (^) and do the same energy estimate as above 
we have the upper bound 

D(^) := lim D(^'") < (|hS"¥). 

A better bound, however, can be obtained for steady, isotropic fiows by a duality argument 
in conjunction with the variational method (see (||12||): 



< C, for a < 



D(f) < V(|hP|2)<| <£p fora>0 (50) 

< i/log (!/£>), for a = 



which agrees with results by other approaches (such as RNG calculation of and Green's 
function method of p7|). 

When the sampling drift is present, the estimates (|38D can be used to compare the 
transport effects of the sampling drift and the subcritical wave numbers. For a + f3 > 1, we 
have from (|38D the bound for the time scale of the fluctuation of particle motion 

e-yD^'\e) » £-V^=("+^-i) = e-2/3/("+2/3-i)^ 

which, in Regime II, is the time scale of observation as determined from the sampling drift 
alone (cf. (P^). Therefore, the transport in Regime II is dominated by the sampling drift. 
For a + 2/3 > 2, (pO]) implies the bound for the time scale associated with the fluctuation 
of particle motion 

which, in Regime III, is the time scale of observation as determined from the sampling drift 
alone (cf. (p4D). Thus, again, the transport in Regime III is dominated by the sampling 
drift. 
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3.3.2 Case 2: Subcritical cutoff 



When the samphng drift is absent, i.e., 7 < 7c, or neghgible, instead of (^8]), we have 

nff < C, ior a + p <l ora<0 



D\f < { 



{ ^£27(1-0-/3)^ for a + /3> 1 
Clog(l/e^), fora + /3 = l 
Ce-T", for a > 



(51) 



^ CJlog(l/£T), for a = 



Vi, j, for some constant C > 0. Note again that e '''" is a better bound than e^'''^^ " ^^ for 
a + 2/?>2. 

The estimates ( |5TD are derived by the same energy estimate as before. In this case, Q = 
and the velocity field consists entirely of the subcritical wave numbers. As a result, <^ in 
(§|) becomes < in (§T|). 

When the sampling drift is absent, the estimates (|5lD yield a lower bound for the scaling 
exponent. For a + /5 > 1,7 < 7c the fluctuation of particle motion is 1/e and is much less 
than 

Thus, we have 

g > l+7-7(a + /5), ioi a + 13 > I. (52) 

The bound ( p^ is sharp when temporal fluctuations of the velocity flelds are the dominant 
mechanism for transport as in Regime V. 

For a + 2/9 > 2, 7 < 7c, (§Tp implies 



^-1 < V£-7«e-2g = ^-g-7«/2_ 

Thus, we have 

g > 1 - 7a/2, for a + 2/3 > 2, (53) 

which turns out to be sharp in Regime VI. 

In the case of subcritical infrared cutoffs, it is often useful to know if the wave numbers 
|k| ~ 5 dominate the transport or not. For this purpose, we modify the previous variational 
method to estimate the transport effect of the wave numbers much larger than the subcritical 
infrared cutoff. We replace c^ by the velocity field Ue consisting entirely of wave numbers 
|k| ~ 5 = e'^ and V^^^ by the velocity field consisting of the wave numbers |k| ^ e'^ . After 
other corresponding modifications are made, the variational method and the subsequent 
energy estimate after substitution of the trivial trial function work the same way. We get 
the upper bound for the contribution, denoted by D*^"^) = [-D>J ], of wave numbers |k| ^ e^ 

to the total eddy diffusivity -D^j • 

A-f < C, for a + /3 < 1 or a < 



w « 



^27(l-a-/3)^ 

log(l/£^),' 

^-70 



l0g(l/£^). 



for a + /3 > 1 
for a + /3 = 1 
for a > 

for a = 



(54) 



Vi, j, for some constant C > 0. For specific application of these bounds, see discussions 
for Regimes V and VI. 

4 Phase diagram for a < 1 

4.1 Regime I: diffusive limit 

First of all, as discussed in Section 3.1, the sampling drift is negligible in this regime: In 
(H), if a + /3 < 1 and a + 2/3 < 2, then 2(3/{a + 2/? - 2) > 2; if a < and a + 2/3 > 2, 
then 2 — a > 2. In either case, the transport would be dominated by the fluctuating wave 
numbers. 

When a + ^ < 1, (|3|) imphes 

< lim inf D^"^ (e) < lim sup D^^^ (e) < cx). (55) 



*o 



*o 



As we will see below that D^^\e) should not vanish in the limit, so the scaling is diffusive 
g = 1 and the scaling limit should be a Brownian motion {H = 1/2). The limit D*{e) = 
liuig^Q D^^\e) , if exists, is the (scale independent) eddy diffusivity. Similarly, for a < 0, the 
variational bound (pO|) implies the diffusive scaling limit. 

The eddy diffusivity probably does not vanish in the limit for the following reason. From 



the turbulent diffusion theorem for mixing flows, proved in |]15|, we know that, for /3 = 0, a < 
1, the scaling is diffusive (g = 1) and the limit is a Brownian motion. As (3 increases, the 
velocity correlation in time increases and so should the rate of transport. But the upper 
bound ( p5D for the eddy diffusivity tells us that it can not enhance transport to the extent 
of changing the scaling limit as long as a + /3 < 1 (This scenerio has been rigorously justified 
in the region a < 0, /3 < 1 in a different turbulent diffusion theorem for non-mixing flows. 



proved in [[15 



For a < 0, the (ordinary) vector potentials for the flows are time-stationary, space- 



homogeneous and have finite moments. Then the diffusion limit theorem of [|T4| holds for 
such flows if molecular diffusion is present (i.e. q = 1,H = 1/2 ii k > 0). And the 
effective diffusivity can be determined from a pair of variational principles ([jl9[)) one of 
which is (IsBf ). This is manifest in the existence of (ordinary) vector potentials when a < 0. 
As shown in [|1^], for steady flows, the existence of space- homogeneous (ordinary) vector 



potentials is the sharp condition for a diffusive scaling limit with molecular diffusion. As 
time dependence of velocity becomes important with (3 < 1, the phase boundary defined 
hy a + P = l,P < 1 points to the fact that the existence of space-homogeneous, fractional 
vector potentials becomes the criterion for the diffusive scaling limit. 

For a + P > l,a < (thus, /? > 1 and a + 2/3 > 2), both the sampling drift and high 
wave numbers are negligible. The effect of molecular diffusion may not be negligible since 
the Laplacian in the Lagrangian generator (|T7|) dominates over A in low wave numbers. 
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4.2 Regime II: space-freezing property 

As we have seen from the analysis of the samphng drift and the apphcations of variational 
bounds, for 

7 > 7c, a + (3>l, a + 2/3<2, a < 1, 

the sampling drift dominates the transport, and, therefore, 

a + 2p — 1 

by (^) (and the discussion afterward). Moreover, since the sampling drift is asymptotically 
uniform in space, the displacement can be approximated asymptotically by 

x%0)+e V(x%0),s)ds. (57) 

Jo 

Eq. ( ^71) is called the space-freezing approximation, in which the space dependence of the 
Lagrangian velocity is suppressed. Eq. (0) defines a Gaussian process with stationary 
increments. It is easy to check that (^) converges to a fractional Brownian motion ^nit) 
by computing its covariance tensor 



2H 



with the Hurst exponent 



{B„(()»B„(()> = C( 



1 a+^_ 1 

2 2/3 '^ ^' 2 



and the coefficient 



For a different and rigorous approach to the fractional Brownian motion limit, see Ref- 
erence 
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4.3 Regime V: subcritical cutoff 

If the cutoff is supercritical, 6 -C k^, the sampling drift is effectively intact, so the frozen 
path approximation (|57D holds along with the FBM limit with q given by (^61). 

If the cutoff is subcritical, S 3> kc, the sampling drift is absent and the transport is 
determined by the fluctuating velocity field. We further decompose the subcritical wave 
numbers into those |k| ~ 5 = e"^ and those much larger. We have made the estimate (([5^) 
for the contribution of the latter to the eddy diffusivity. 

By a simple spectral calculation, the velocity field \Je{'^, t) consisting of the wave numbers 
|k| ~ 5 can be approximated by 

5^-"U(5x, ^^/^t) (58) 
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where U has the energy spectrum (^ for |k| G [1,C] with a sufficiently large constant C. 
Substituting (|^) into the equation of motion we have 

dx' = e^(i-")-2g+iu(a;=(t)/ei-7^ t/e^''-^^^)dt. (59) 

The effect of the wave numbers |k| ^ 5 is like adding a turbulent diffusivity to eq. (|5^), i.e. 



dx' = e^(^-")-2y+^U(a;"(t)/£^-^, t/e^''-^^^)dt + e^~''^2'b^^)d'B{t) (60) 

where B(t) is the standard Browian motion and D^^-* is the portion of the eddy diffusivity 
coming from the wave numbers |k| ^ 5 (cf. the discussion preceding (plD). 

Since U is a mixing flow, we expect the limit of (^) to be a Brownian motion. We also 
expect the time variable in (|59D to dominate, so we equate 57(i-a)-2g+i _ ^213-2/37 ^^^ arrive 
at the expression 

g = l + 7-7(a + /3) (61) 



as one would expect from a generalized 'diffusive' scaling of (|59|). We check that the space 
variable in (|59|) is indeed relatively slow in the sense 

(l-7)/(g-/57)<l, (62) 

for a + 2/? < 2. With (|TD, eq. (||) becomes 



dx' = r/-iU(a;/r/f ^^)/(«-^^),t/r/2)t/t + e^-''\/2t>(^)dB{t), with r?, = e""^^. (63) 



The bound (|5|) and §^) imply that e^~W2t>(^) < ^i-g+7-7(«+/3) = 1 as e tends to zero. 
Subcriticality, 7 < l/(a + 2(3 — 1), implies r/^ -^ 0. Eq. ( |63D satisfies the conditions of 
the diffusion limit theorem of ^^ and |^ for mixing flows with a generalized 'diffusive' 
scaling ([63|)-(|6^). The limit is a Brownian motion with the diffusion coefficients given by the 
Eulerian Taylor-Kubo formula 

/ (U(0,t)®U(0,0))rft. (64) 

Jo 

The condition a + (3 > 1 implies g < 1. 

The scaling law with H = 1/2, q, given by (0), and the eddy diffusivity given by the 
Eulerian Taylor-Kubo formula (^) holds in the other part of Regime V (a > 1) as well (see 
Section 5). 

4.4 Regime II': critical cutoff 

When the infrared cutoff is critical i.e. 7 = 7c = {a + 2(3 — 1)^^ the critical wave numbers 
are still present. Therefore the scaling exponent is given by (^). 

Rescaling the velocity field XJ^ consisting entirely of wave numbers |k| ~ e'^'' but |k| < e^" 
by (p8|), we have, instead of (|60D , the equation 

dx' = IJ {x'{t)/e^~^%t)dt + £^-V2D(=)(iB(t) 
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where U has the energy spectrum (|28|) supported in [1, oo) and D^^^ is the eddy diffusivity. 
Since 7c > 1 and g < 1, we have in the hmit e —>■ 

dZ{t) = V{0,t)dt. (65) 

Because the energy spectrum of U does not have small wave numbers, the process Z is 
probably not self-similar. Thus, the Hurst exponent is not well-defined. However, the long 
time asymptotics of (|65|) is a Brownian motion, due to the mixing property of U, so we may 
associate the asymptotic Hurst exponent 1/2 to the process defined by (|65D. 

4.5 Regime III: time- freezing property 

In this regime, the sampling drift is a time dependent, mean zero random variable whose 
second moment is of order 

' " |1-2q: Till 2-2a 






The mixing time for the sampling drift is no less than 1/e^^. On this time scale the sampling 
drift transports particles over the distances no less than e^~'^e~'^^ = gi-^-^/? -y^^j^j^j^ jg no less 
than the spatial observation scale 1/e if a + 2/5 > 2. Thus, the sampling drift appears steady 
on the observation scale. The time to exit a ball of radius 1/e, based on the sampling drift 
alone is e~^e"~^ = e"""^. 

From ( pO] ) it follows that the time to exit a ball of radius 1/e, with its center moving 
by the sampling drift, is ^ e"~^, as a result of the spatially fluctuating, subcritical wave 
numbers (|k| ^ e). Thus, by (P^, the effect of the subcritical wave numbers is dominated 
by that of the sampling drift. Therefore, by (0), we have 

q=l- a/2. (66) 

The scaling exponent ( p^ ) is superdiffusive for a > 0. Note that 

q - P/{a + 2/3 - 1) = (a/2 + /?-!)(!- a)/{a + 2/5 - 1) < 

for 1 < a + 2/? < 2, a < 1 and, thus, for the same a, the rate of transport in Regime II is 
smaller (since faster decorrelation in time tends to slow down the transport). 

Since the transport is dominated by the effectively steady sampling drift, we may consider 
the velocity field XJg consisting entirely of the wave numbers |k| ^ kc = e and freeze the time 
variable in the resulting velocity field as e — *> 0. Rescaling as in (|58|) , we have the equation 



dx%t) = e^-''-^'^lJ{x%t), 0)dt + e^~''V2i^dw{t). (67) 

Eq. (^Tj) has a limit if and only if (|66|) holds, which also make the diffusion term vanish in 



the limit for a > 0. Formally, the limit process Z satisfies 

dZ{t)=V{Z{t),0)dt (68) 

where U has the energy spectrum (^ supported in |k| G (0, oo). The supercritical cutoff is 
now removed by the rescaling (|58D. This indicates the limit Z is self-similar. 
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It should be noted that the velocity field U is a generalized function because of the 
ultraviolet divergence in the energy spectrum of U. Thus, eq. (^) is not well-defined in the 
ordinary sense. Study of transport in generalized velocity fields is interesting by itself, but 
we do not pursue it here. Our only purpose is to use the energy spectrum of the velocity 
field as an indicator of the self-similarity of the limit process and to show how the space and 
the time dependence of the velocity field enter the equation. The same remark applies to 
the same situation in the sequel as well as (|65D and will not be repeated. 

We now identify the Hurst exponent of Z. From (H), we have the asymptotics for the 
covariance of the successive increments on the time scale t ~ e^'^'^ 

{{X{2t) - X{t)) ■ {X{t) - X{0))) ~ t^H _ ^-4gif_ ^gg) 



On the other hand, by (|2^ , the covariance of the successive increments on the time scale 
t ~ e~'^'' is of order (/c^'^e"^"?)^. Equating it with (p^D, we have H = l/{2q). We hypothesize 
that the limit be a fractional Brownian motion. 

Molecular diffusion probably has no significant effect on the scaling law, even though the 
presence of molecular diffusion is needed in the variational principle (0). As remarked after 
(|T7|) in Section 2, the molecular diffusion is negligible for (3 < 1. For /3 > 1, as larger (3 
gives rise to longer correlation times, the scaling law should have equal or smaller scaling 
exponent q. However, since the time-freezing property has already set in for /5 < 1 and 
resulted in a scaling exponent independent of P, the absence of molecular diffusion would 
not have changed the scaling exponent for /3 > 1. 

4.6 Regime VI: subcritical cutoff 

The results of the previous section hold for any supercritical cutoff (7 > 7c = 1) as the 
sampling drift is essentially intact and dominates the transport. 

We separate the wave numbers |k| ~ 5 from |k| ^ 5 and rescale the equation as in 
Section 4.3 to obtain (|60D . As before, we expect the limit to be a Brownian motion. In this 
case, however, we expect the space variable in the velocity field to dominate the transport. 
So we equate £'>'(^~")~29+i = ^7-1 g^^^d obtain the scaling exponent 

q = l- 7a/2. (70) 

Rewriting (^) with r]^ = e^~'^ we have 



dx%t) = ri;^Vix'{t)/rie, t/r^f «-^^)/(i-^V^ + e^'V 2D(^)rfB(t) + e^~^V2i^dw{t). (71) 

The bound (||) and (^ imply that e^^'^ViD^ < ^1-9-70/2 ^ i^ gjj^pg g < 1 the factor in 
front of the molecular diffusion also tends to zero as e ^ 0. The time variable is relatively 
slow in the sense 

as a + 2^ > 2. 

The velocity field U has a fixed infrared cutoff and, consequently, gives rise to a space- 



homogeneous vector potential, so, as suggested by the diffusion limit theorem of |jTj], the 
limit should be a Brownian motion. 
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The diffusion limit theorem of |14], however, does not apply directly because it was proved 



with a non-vanishing molecular diffusion. Generalizing the theorem of |0 to the situation 
with a vanishing molecular diffusion such as eq. (|7lD remains a challenging problem in 
turbulent transport. 

4.7 Regime III': critical cutoff 

When the cutoff is critical, i.e. 7 = 1, the critical wave numbers are still present and 
dominates the transport. Therefore, the scaling exponent is given by (|66D. 



The difference is that eq. ([67| ) now has the limit satisfying eq. ( |68D with the energy 
spectrum of U supported in |k| G [1, 00) rather than (0, 00). As a consequence, the limit is 
not self-similar. 

4.8 Phase boundary 

The transport for the phase boundary a + 2(3 = 2,0 < a < 1 contrasts interestingly to 
regime on either side of the boundary. 

When the cutoff is supercritical, the sampling drift is present and dominates the trans- 
port. Thus, we expect q = (3 from (0). Indeed, with that and 7c = 1, we have as before the 
asymptotic equation 

dx' = V{x'{t),t)dt + e^-^Vmi^dBit) 



where U has the energy spectrum ( P^ ) supported in |k| G (0, cxd). In the limit, the diffusion 
term vanishes (cf. (|38|) ) and the limit Z satisfies the equation 

dZ{t) = V{Z{t),t)dt (72) 

whose solution is expected to be a fractional Brownian motion. 

When the cutoff is subcritical, we expect g = 1 — 7 -|- 7/? from (^). With that and 
7]^ = e^~'^, we have the asymptotic equation 



dx' = T]-^\J{x/7],,t/rjl)dt + £^"V2D(^)rfB(t), (73) 

where the energy spectrum of U is supported in |k| G [1, 00). The diffusion term dies out in 
the limit as before (cf. O)) and the limit Z is a Brownian motion due to the spectral gap 
in U by a turbulent diffusion theorem of |]T3|. 



Contrary to Regimes II, IF, V, III, III' and VI, both the time and the space dependence 
of the velocity field in ([7^) and (73) affect the transport. 



5 Phase diagram for a > 1 

If the infrared cutoff threshold 5 is fixed as e tends to zero then the velocity is mixing in 
time and, by the turbulent diffusion theorem of |T^ , the scaling limit is a Brownian motion 
{H = l/2). 

Anomalous scaling limits arise when 5 is coupled to the spatial observation scale: 5 = 
£'•',7 > 0. The exponent 7 characterizes the relation between the spatial observation scale 
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Xje and the energy containing scale 1/5. And, in this case, super diffusive scahng results 
from divergent mean kinetic energy as the infrared cutoff is removed; in particular, when, 
the cutoff is supercritical, 7 > max {(a + 2(3 — 1)~^, 1}, (Regime IV), the energy containing 
scale is larger than the spatial observation scale, and it results in a super-ballistic scaling 
q < 1/2 and a regular limit {H = 1). 

Since the transport is dominated by wave numbers |k| ~ 6, it is natural to rescale the 
velocity V as 

V(x,t) = (5^-"U((5x,(52^t) (74) 

where the velocity field U has the energy spectrum (|^ supported in |k| G [1, 00). Contrary 
to U occurring in the case of a < 1, the velocity field U in the case of a > 1 is an ordinary 
function, since there is no ultraviolet or infrared divergence, and temporally mixing due to 
the spectral gap in U. 

In terms of U, the equation of motion becomes 

Now that the infrared cutoff of U is 1, the limit is expected to be a Brownian motion as 
long as either space or time variable is fast. Depending on the parameters, two types of 
diffusion limit theorems in the literature are pertinent to this limit: one is based on velocity 
decorrelation in time (|3^,[2^|, |^, |SC]) and the other on velocity decorrelation in space 



(H)- 

Eliminating the infrared divergence by rescaling is also the approach of Avellaneda and 
Majda ||^, in which the case of a critical cutoff 7 = 1 was considered in the region (3 < 
1/2, < a < 2. Here we adopt the same idea of rescaling and generalize their results by 
using new limit theorems which were not available to them. 

5.1 7 > 1: Regimes V, IV and IV' 

For 7 > 1, eq. (^) does not have fast space variables. To have a nontrivial limit, we must 
have 2g — l + 7(a — 1) > 0: For 7 > 1 space variable is not fast in ([751). To have a nontrivial 
limit, we must have 2g — 1 + 7(a — 1) > or 2g — 1 + 7(a — 1) = 0. The former case gives 
rise to Regimes V wheres the latter gives rise to Regimes IV or IV'. 

5.1.1 Regime V 

When 

2g-l + 7(a-l) >0, (76) 



, U in ([75| ) has a large multiplier, so a nontrivial scaling limit requires rapid time relaxation, 
i.e., 

q > 7/3. (77) 

By choosing a generalized 'diffusive' scaling for eq. ([75|), i.e., ((^""^e^'^^^)^^ = S'^ /e'^ or 

g = l + 7-7(« + /5), (78) 
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([75|) becomes 

dx'{t) = r];^V{4^-^^/^''-'^^^x%t),t/r]'^^)dt (79) 

with ?7£ = e'^~^'^. Eq. ( [7UD has the form of the classical diffusion limit theorem (|3^, |0, P^ ) 
(Moreover, the velocity field U is smooth and satisfies the mixing condition of Rosenblatt 
( |[40|| ) even for /3 > since U has no small k components) Thus, the process x'^{t) converges to 
a Brownian motion {H = 1/2) with diffusion coefficients given by the Eulerian Taylor-Kubo 
formula (|6^). 

However, there is one constraint to be considered: (^) must be consistent with (^), i.e., 
?7e must tend to zero with e. This means 7 < l/(a + 2j3 — 1), a subcritical cutoff. 

5.1.2 Regime IV: smooth motion 

If the cutoff is supercritical, as discussed in Section 3.1, the transport in this regime is 
dominated by the sampling drift that is, in turn, dominated by the wave numbers near by 
the infrared cutoff. Time as well as space dependence of the velocity field are irrelevant. 
Because both space and time variables are slow in the velocity field, nontrivial scaling limit 
holds only if 5"~^e^''~^ = 1. Thus we have 

g = (l + 7)/2-7a/2. (80) 

Consistency, < g < 7/5, then implies that 7 > l/{a + 2/3 — 1) and 

a < 1 + 1/7. 

The limit process Z(t) is advected by a constant drift 

dZ{t) = V{0,0)dt (81) 

and is regular, or smooth {H = 1). 

We note that the limit ( [8T| ) is independent of the initial condition x'^{0), is self-similar 



and has a well defined Hurst exponent. 

5.1.3 Regime IV': critical cutoff 

When 2g — 1 + 7(a — 1) = 0, or, equivalently, 

g=(l + 7)/2-7«/2, (82) 

and 

q = jP (83) 

a nontrivial limit results. Combining ( ^2|) and (|83|) , we have a + 2(3 = 1 + I/7 or 

7 = l/(a + 2/5-1) (84) 

which defines a critical cutoff for a + 2(3 < 2. The limit Z(t) satisfying 

dZ{t) = V{0,t)dt (85) 

26 



is a smooth, Gaussian process. The Hurst exponent is not strictly well-defined for Z(t) due 
to lack of self-similarity. On large time scales, however, (^) has a Brownian motion limit 
as U is temporally mixing and, thus, an asymptotic Hurst exponent H = 1/2. This case is 
similar to Regime V. 

In particular, when 7 = 1 = l/(a + 2/5 — 1) and q = /3, eq. (|75|) is independent of e, i.e. 
x%t) = Z(t) with 

dZ{t) = V{Z{t),t)dt. (86) 

The Hurst exponent is not strictly well-defined for Z(t) due to lack of self-similarity. But, 
as the velocity field is temporally mixing, the long-time limit of Z(t) is a Brownian motion, 
by the turbulent diffusion theorem of ||T^. So if = 1/2 is the asymptotic Hurst exponent. 
Another critical cutoff is 7 = 1 for a + 2/3 > 2. The equation of motion (^) becomes 

dx%t) = e2-°-2'?u(a;"(t),t/e25-2/3)^i 
which has a nontrivial limit when q = (3 = 1 — a/2. The limiting process Z(t) satisfies 

dZ{t) = V{Z{t),0)dt (87) 

which is regular for finite times, but is not self-similar. Thus the Hurst exponent is not well 
defined. It is not clear whether an asymptotic Hurst exponent is well-defined either, since 
we do not know if, without molecular diffusion, motion in three dimensional, steady flows 



like (|87|) can be homogenized or not. Bounded and unbounded streamlines may co-exist in 



steady flows and, if so, the resulting limit would depend on initial conditions. 

5.2 7 < 1: Regimes V and VI 

With 7 < 1, fast space variables now enter the picture. There are two regimes depending on 
whether time relaxation dominates over space decorrelation. 

5.2.1 Regime V: dominant time relaxation 

This regime occurs when 6^ /e'^ ^ 6/5 or, equivalently, 

g > 1 - 7 + 7/3. (88) 

Then, by choosing the 'diffusive' scaling, (^""^e^g-i^-i = S^/e*^, i.e., 

g = l + 7-7(a + /3), (89) 



eq. (|7^) can be rewritten as ([79|), except that the space variable is also fast, albeit not fast 
enough to have an impact in the diffusive scaling. In this case, the generalized limit theorems 
proved in [^,|3^ are applicable and they extend the validity of the Taylor-Kubo formula to 
our situation. (Like the classical diffusion limit theorem, the generalized limit theorems also 
require the mixing condition and regularity on the velocity U, both of which are satisfled 
here.) The limit is a Brownian motion {H = 1/2) with diffusion coefficients given by the 
Eulerian Taylor-Kubo formula (|6^). 



Condition ( pH]) requires that 

a + 2/5 < 2. (90) 
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5.2.2 Regime VI: dominant space decorrelation 

For 

g< 1-7 + 7/3 (91) 

velocity dependence on space now dominates over the dependence on time in the diffusive 
scaling of eq. (|75[) . By choosing the scaling (^"~^e^^~^)~^ = 5 /e or, equivalently, 

g = 1 - 7a/2, (92) 

Eq. (|75[) is rewritten as 

dx'{t) = r];^Vix%t)/rie, t/r^f '?-^^)/(^-^))rft, with r/, = e'-^. (93) 



The limit of (|93|) should be a Brownian motion as U gives rise to a space-homoegeous vector 
potential (cf. the discussion in Section 4.5). 

Consistency ((0) and q > 0) requires that a + 2/3 > 2, 7 < 2/a. 

5.2.3 Phase boundary 



On the phase boundary a + 2(3 = 2, eq. (p3D becomes 



dx%t) = 7];^lJ{x%t)/7],, t/r]f)dt (94) 



with a temporally mixing flow U. Space and time correlations play comparable roles in 
From the turbulent diffusion theorem for mixing flows ([^), it follows that the solution 
x^{t) has a Brownian motion limit {H = 1/2). 

6 Conclusions 

The supercritical and subcritical diagrams (Fig.l and Fig. 2) are divided by the line, a+/3 = 1, 
and the line, a + 2/3 = 2, and/or the vertical lines a = 0,1,1 + I/7, 2/7. First of all, 
a + (3<lora<0 defines a cutoff independent diffusive regime, in which the sampling drift 
is negligible. Outside of the diffusive regime, the line, a + /3 = 1, is the cross-over between 
short-ranged and long-ranged velocity correlations; the latter manifests in the fact that the 
sampling drift dominates the transport and subcritical wave numbers are negligible. The 
line, a + 2/3 = 2, is the cross-over between velocity dependence on space and on time; in the 
region above the line, velocity dependence on time is negligible, whereas in the region below 
the line velocity dependence on space is negligible. 

Fig. 3 and Fig. 4 are cross sections of the full three-dimensional phase diagram at 7 = 
const. > 1 and 7 = 1, respectively. In Fig.3, the cutoff is supercritical for a + 2/3 > 1 + I/7 
and subcritical for a + 2/3 < 1 + I/7. In Fig. 4, the cutoff is subcritical for a + 2/3 < 2 and 
critical for a + 2/3 > 2. Cutoffs with 7 < 1 are subcritical and, thus, covered in Fig. 2. 

The limit is one of the three kinds: Brownian motion {H = 1/2), persistent fractional 
Brownian motion (1/2 < H < 1) or regular, or smooth, motion {H = 1). The relation 
H = l/(2g) holds for a < 1 with supercritical infrared cutoff but neither for subcritical 
cutoffs nor for a > 1 (In these situations, H < l/{2q), instead). For the critical cutoff 7=1, 
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the Hurst exponent is not well defined. However, an asymptotic Hurst exponent may be 
defined and it is equal to 1/2. The diffusive regime (g = 1, if = 1/2) is most robust in that 
the scaling law is independent of any infrared cutoffs. The fractional Brownian motion limit 
of Regime II and HI is not affected by supercritical cutoffs. All other regimes are cutoff 
dependent. 

In the case of subcritical infrared cutoffs, with the rescaling of the velocity field, the dia- 
gram can be understood by means of three types of diffusion limit theorems in the literature: 
(i) one for which the spatial dependence of velocity is negligible and the effective diffusivity 



is explicitly given by the Eulerian Taylor-Kubo formula (|^5[,[gJ],[^,|3^), (ii) another for 
which the temporal dependence of velocity is negligible, but molecular diffusion is assumed 
to be present, and the effective diffusivity is implicitly given by a pair of variational principles 
([pTSf, [0], |jT9|) or a Stieltjes integral formula ([§]), and (iii) the other for which space and 



time dependence of velocity play comparable roles and is referred to as turbulent diffusion 
theorems (Two such theorems are proved in |1^). As in (ii), the turbulent eddy diffusivity 



can be written as a variational principle similar to (0). After rescaling. Type (i) limit theo- 
rems apply to Regime V; Type (ii) limit theorems apply to Regime VI and part of Regime I; 
Type (iii) limit theorems apply to Regime I or on the phase boundary a + 2/3 = 2, < a < 2. 
All three types of limit theorems are insensitive to the dimension. This may explain certain 



similarity between the phase diagram for shear-layer flows of [|l| or |^ and the diagram for 
isotropic flows with a subcritical cutoff. 

Similarly, there should be three types of fractional-Brownian-motion limit theorems: one 
completely determined by Eulerian time decorrelation (Regime II), another determined by 
the space decorrelation (Regime III) and the other determined by the dependence of velocity 
on both space and time (The phase boundary a + 2/5 = 2, < a < 1). 

In the case of supercritical infrared cutoffs new phenomena emerge: dominant sampling 
drift, fractional Brownian motion limits, critical infrared cutoff and related cutoff depen- 
dent effects. Although these phenomena are introduced and analyzed for motion in three- 
dimensional, isotropic flows, they also arise in two-dimensional flows or anisotropic flows 
such as random shear-layer flows. A difference lies in the role of molecular diffusion which 
is much more prominent for anisotropic or two-dimensional flows (see discussion in the next 
section). New variational principles for the cutoff dependent eddy diffusivity are formulated 
and used to obtain general bounds for the eddy diffusivity. 

Contrary to subcritical and supercritical cutoffs, regimes (IF, HI' and IV') with criti- 
cal cutoffs produce limits that are not self-similar and do not possess a well-defined Hurst 
exponent. 

Scaling limits of turbulent transport in flows with a nonzero mean drift have different 
phase diagrams (see |2^, |^, 0) and will be reported in a forthcoming paper. 



6.1 Role of molecular diffusion 

Molecular diffusion has at least two roles: (i) to eliminate possible dependence of scaling limit 
on the initial point, as particles may be trapped by closed or bounded streamlines, so that 
the process may be homogenized; (ii) to reduce dynamic velocity correlation in time and, 
thus, change the scaling law to one with larger scaling exponent q. The first role of molecular 
diffusion is prominent for transport in steady flows; without it, localized streamlines would 
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prevent homogenization from happening (see [18|). In this connection, molecular diffusion 
also helps blend the effects of streamlines of different scaling behaviors. In the present work, 
we always assume homogenization and we focus on the second role of molecular diffusion. 

Molecular diffusion is negligible in Regimes II, III, IV, where the sampling drift dominates 
the transport as well as Regimes V and part of I (i.e., a + (3 < 1 in the limit of high Peclet 
number), where velocity decorrelation in time is significant. But its effect is not so clear in 
Regimes VI and the other part of Regime I (i.e., a + P > l,a < 0), where the spatially 
fluctuating wave numbers dominate. 

In this regard, when high wave numbers are negligible, as for a+jS > 1, one can go further 
by comparing the term representing molecular diffusion, kA, and the term representing the 
flow. A, in (p!7|), and see that, for /3 < 1, the effect of molecular diffusion should not affect the 
scaling law. As the velocity dependence on space dominates over that on time, the scaling 
law is independent of (3. In view of the discussion in Section 4.4 on effect of molecular 
diffusion, we expect the scaling law for (3 > 1 can be extrapolated from that for j3 < 1 to 
conclude the scaling law of Regime VI is independent of molecular diffusion as well as (3. 

As a on the phase boundary (a = 0) is bigger than that in the region (a + /? > 1, a < 0), 
for given (3, the scaling exponent q in the region, with the absence of molecular diffusion, 
should not be less than that of the phase boundary, which is 1. Therefore, g = 1 in this 
region even without molecular diffusion. 

7 Appendix: variational principle for eddy difFusivity 

To derive the variational principle (^) we consider a pair of period cell problems for an 
arbitrary constant unit vector e 

^fen)^J.-) + (c^ + V(",n)).VxJ'"^+V(^'")-e = 0, inf2(") (95) 

^(^'"V-'"^ - (c. + V(^'")) • VxL"'"^ - V(^'") • e = 0, inO(") (96) 

where both x+ ^^^ X- satisfy the periodic boundary condition. Note that (^) is simply 
the adjoint of ( pS] ) as V^*^'"^ is divergence free. 

Adding and subtracting ( p5D and (p^) we obtain 

^ie,n)^(e,ny ^ ^^^ ^ y(e,n)^ . y^(.,n) ^ y(e,n) . e = Q (98) 

where 

X^^'") = ^(xS'"^ + X^''% X^^'"^' = ^(xi^'"^ - X^^'"^ (99) 

First we establish some useful identities for 

D(^'")(e) = D(^'")e ■ e = -(^(^'"^x+'^^xS' "V- 

Proposition 1 

D(^'")(e) = (xj'")v(^'") ■ e). = -(xL^'")v(^'") ■ e). = -(^(^'")x^''"^X^''"^)n. (100) 
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The first identity in ( [1001 ) follows from integration by parts after multiplication of eq. (p5| ) 
by x+ ■ To verify the second, we make use eqs. (|95|) , (|96D and the divergence free property 
of V*^^'"-* in the following calculation 

Here we have used the identity 



C, + V(^''^)) ■ VxV^'"^J Xr% = =f ■ {{Ce + V(^'"))(xV^'"T)n = 

as a result of the incompressibility of c^ + V^^'"^ and the space- homogeneity of ((ce + 

v(^-"))(xS''^))2)„. 

Thus, in view of (PUJ), the following result is clear. 
Proposition 2 

D(^'")(e) = (;^(^.")'v("'") ■ e)„ = -(^(^'")x^''"^'x^''"^')n - (V("'") ■ Vx^"'"^X^''"^')n- 

Next we derive the variational principle (^Oj). 

Let g be the minimizer of the convex functional in (|30|) and (?' be the periodic solution 
of the equation ^ ^ 

^("■")^' + (c, + V("'")) ■ V^ + V("'") -6 = 0. (101) 

Taking the first variation of the functional in ( pOD at g we have 

- (^(^■")^5^)„ - {A^''-'>g'6g')n = (102) 

where the variation 6g' is related to the variation 6g by 

^("'") V + (c, + V(^'")) ■ V5^ = (103) 

following (|101|) . Substituting (|103|) into (|102|) and integrating by parts we get 

(^^^'"^^5^)n +{{cs + V(^'")) ■ Vg'6g)n = 

for all admissible variations 6g. Thus 

A^''"'^g + (c, + V(^'")) ■ Vg' = 0. (104) 

Since eqs. (0) and ( pSj) (also ( |101| ), (|104| )) are well posed, we conclude that g = x^^'^''\g' = 
■^ie,n) ^p ^Q constants. 

By reversing the above argument, it is easy to see that g = x^^'"'\fl'' = x*''^'"'' with 
^(e,n) ,^(£,n)' gjyg^ }^j ^ggj) g^^e the minimizer of (|30|). 
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Guide to the phase diagrams 

The full diagram is in the three dimensional space of a, /3, 7. To simplify the presentation, 
we choose to portray the full diagram as four planar diagrams: the supercritical diagram, 
the subcritical diagram, the cross sections 7 = const. > 1 and 7 = 1. Note that Regime IF 
does not show up in any of the diagrams. 

• The exponents a and (3 defined in (§)-(0) characterize the space-time correlations. 

• The exponent 7 is related to the infrared cutoff 6 in (^ as 6 = e'^. 

• The critical exponent 7c = max {1, (a + 2/? — 1)^^}. 

• The scaling exponent q and the Hurst exponent H are defined by (H) and (|), respec- 
tively. 

• Regime I: a + /? < 1 or « < 0. 

Regime II: a + /? > 1, a + 2/? < 2, a < 1 with 7 > l/(a + /?-!). 

Regime II': a + (3 > l,a + 2(3 < 2,a < 1 with 7 = (a + 2/3 - 1)"^ 

Regime III: a + 2P > 2,0 < a < 1 with 7 > 1. 

Regime III': a + 2(3 > 2,0 < a < 1 with 7 = 1. 

Regime IV: I < a < 1 + I/7 with 7 > max {1, (a + P - l)^^}- 

Regime IV': 1 < a < 1 + I/7 with 7 = 7c. 

Regime V: a + (3 > l,a + 2(3 < 2 with 7 < l/(a + (3-1). 

Regime VI: a + 2(3 >2,a>0 with 7 < 1. 
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Fig.l. Phase diagram with supercritical cutoffs : y > max{l/(a+2p-l), 1 } 
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Fig. 2. Phase diagram with subcritical cutoffs: y< max{ l/(a+2p-l), 1 } 
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a=l 



a = 1 + 1/y 
Fig. 3. Phase diagram at any cross section y = const. > 1. 
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Fig. 4. Phase diagram at the cross section y = 1. 
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